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Abstract 

By solving exact three-body equations, we study the three-baryon system with charm -|-1. We 
look for possible bound states using baryon-baryon interactions obtained from a chiral constituent 
quark model. The smaller effect of the Ac -H- Sc conversion reverses the order of the (/, J) = (0,1 /2) 
and (/, J) = (0, 3/2) states, rather close on the strange sector. The diminishing of the kinetic energy 
due to the large reduced mass gives rise to a bound state in the {I,J) = (0,3/2) channel. After 
correcting for Coulomb effects the binding energy would be between 140 and 715 keV. 
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Soon after the discovery of baryons possessing net charm it was suggested that there 
should also exist charmed nuclei. The observation of a candidate event that could be in¬ 


terpreted in terms of the decay of a charmed nucleus [l|, fostered conjectures about the 
possible existence of charm analogs of strange hypernuclei j2“4|. Three ambiguous can¬ 
didates of charmed hypernuclei were reported by an emulsion experiment with 250 GeV 
protons j^. This gave rise to several theoretical estimates about the binding energies and 
the potential-well depth of charmed hypernuclei based on one-boson-exchange potentials for 
the charmed baryon-nucleon potential jo-lO]. There were also theoretical estimations of 


the production cross-sections as well as experimental condition requirements for producing 
charmed hypernuclei by means of charmed exchange reactions on nuclei [n|. The experi¬ 
ments for searching charmed hypernuclei are becoming realistic and may be performed in 
coming yea rs at Hall C of JPARC |^, at the FAIR experiment 13|, or at the SuperR 


collider 


14| . All these experimental prospects have reinvigorated twenty years later the 


study of charmed hypernuclei ISHlTj and also more recently theoretical studies of charmed 
dibaryons 18|, ll9|. We show in Fig. [Da diagram to generate A^ baryons by means of an¬ 
tiproton collisions on the deuteron through the intermediate production of charged D mesons 
(m£)± = 1869.61 MeV/c^) feasible at modern factories, that has been proposed to study the 
existence of charmed hypernuclei at JPARC 12|]. A similar reaction, + p ^ A+ + tt’*'. 


is proposed at SuperR to detect charmed super-nuclei 


1 ^. 


We have developed in the past the exact formalism to study three-baryon systems with 
a heavy flavor baryon |20|. It is our purpose in this work to extend our previous study to 
three-baryon systems with a unit of charm by considering the region of the AcNN bound 
states, where a charmed hypertriton might exist, and calculating for the hrst time the A^d 
and Scd scattering lengths. We will simultaneously study all AcNN and HcNN states with 


vP 

N. 

D“ 

Tl 

FIG. 1: Diagram to generate by means of antiproton collisions on the deuteron. 
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J = 1/2, 3/2 and I = 0,1,2. 

Let us start by summarizing the description of the three-baryon system with a charmed 
baryon. We transform the Faddeev equations from being integral equations in two continuous 
variables into integral equations in just one continuous variable by expanding the two-body 
t—matrices in terms of Legendre polynomials. Pm 2l| , 


ksidPuP'iie) = J2Pn{Xi)Tl^^,.i.{e)Pr{x'i), 


( 1 ) 


where x* = x' = Pi and p' are the initial and final relative momenta of the pair 
jk, and 6 is a scale parameter of which the results do not depend on. 

If we identify particle 1 with the charmed baryon and particles 2 and 3 with the two 
nucleons, the integral equations for pd scattering at threshold, with /? = or A^, are in the 
case of pure S'—wave configurations, 

P+(T1+(T3 — S2+T1+T3 —*2 


ms^is 


rTinS 2 l 2 f ^ \ _ JDnS 2 l 2 (^ 

/•OO 

+ 2J2 dq,Alifr'(q,,quE)Ari:sT‘"’(quq,\E) 
Vo 


rsiii 


T^riQs), (2) 


where ai (ri)and (ts) stand for the spin (isospin) of the charmed baryon and the nucleon 
respectively, while Sj and F are the spin and isospin of the pair jk. T 2 ^gj^j^{q 2 ) is a two- 
component vector, 

(T2-SI% 0i(l2)\ 

rs’y*) = dnlF . 


\T?i%,q(q2) 

while the kernel of Eq. ([2]) is a 2 x 2 matrix defined by 

/ qs; E) qA E) \ 

V gs; E) A-™(.^2, gs; E) J ’ 

/ ^ 31 ;S/;ScJV(Sc) (q2,qi;E) A-^^;^-^^^^Jq2,quE) 

^ ^31;S/;a1v(Sc) (q 2 ,quE) 

/^13;S/;Af|;!(gi) g3j-^) 0 \ 

V 0 Arsl^^-l^(quqs;E)J ’ 

where 


(3) 


A-gr^^Hq2,q3;E) = 
Asi$r^^Hq2,qi;E) = 

^.™«3(gi^g3.^) ^ 


( 4 ) 


gi 


AT3$}7riq2,q3-,E) = ^ - g2V2^2)f 


nl 


X 


Pr{x'2)Pm{x3) 


/ dcos9 -— ^ , 

1-1 E + AE6pA^ - P3/2/i3 - g|/2z/3 -i- le 


, o, /3 Eg, Ag, 
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9? 


r ^ 

Pr{x'.;^Pm{X\) 


X / (icos 6 ' —-—„ ^ , 

7-1 E + Ai?(5/3A^ - pf/2/ii - g^/2z/i + ze 


, Q^, /3 Sc, Ac, 


21?3”s™”(®. ®; -B) = ASS?" E (B + ABi^A. - 9 iV2i/i) 

T ^ 

Pr{x'^Pm{x-i) 


X 


(icos 6 * ■ 


-; /? = Sc,Ac, (5) 


J-1 E +-P3/2/i3-g3/2z^3 + *e’ 

with the isospin and mass of particle 1 (the charmed baryon) being determined by the 
subindex /9. /ij and z/, are the usual reduced masses and the subindex aN[(3) indicates a 
transition aN fiN with a nucleon as spectator followed by a NN —)• NN transition 
with j3 as spectator. '7‘^52i2;«/3(^) coefficients of the expansion in terms of Legendre 

polynomials of the charmed baryon-nucleon f—matrix t 2 ;s 2 * 2 ;a/ 3 (P 2 ,P 2 ) ^^e transition 

aN —)■ f3N^ i.e., 


fpi = 


2?7, + 1 2r + 1 P 


1-1 


dXi / dx[Pn{Xi)ti.siii;al3{Pi,Pi,e)Pr{x'i} . 


1-1 


( 6 ) 


The energy shift AE, which is usually taken as — Ms, will be chosen instead such that 
at the 13d threshold the momentum of the ad system has the correct value, i.e., 

[(m^ + mdY - {rria + rndf] [{mp + rUdY - (m„ - mdY] 


A£/ = 

^Padimp + rUdY 

where pad is the ad reduced mass. 

The inhomogeneous term of Eq. ([2]), i? 2 ; 5 p/ 3 ('? 2 ), is a two-component vector 


(7) 


B2-JYi3{Q2) — 




where 


7?7;5Pa/l(g2) = - ql/2l^2)Pr{x'Y. 


( 8 ) 


(9) 


with Si = 1 and A = 0 are the spin-isospin transition coefficients corresponding to 
a charmed baryon-deuteron initial state (see Eq. (30), of Ref. |2l|), 0d(?2) is the deuteron 
wave function, is the energy of the /3d threshold, Pr{x 2 ) is a Legendre polynomial of 
order r, and 

—q2 — b 

m 3 
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Finally, after solving the inhomogeneous set of equations ([2]), the 13d scattering length is 
given by 

( 11 ) 

with 

fOO 

Tw = 2 E fclS? / qldqM']^)Pn{x'PnS;UaP- (12) 

123212 

In the case of the Tj^NN system, even for energies below the Serf threshold, one encounters 
the three-body singularities of the A^iViV system so that to solve the integral equations ()2|) 
one has to use the contour rotation method where the momenta are rotated into the complex 
plane qi —)■ qie~'^3 since as pointed out in Ref. [2l| the results do not depend on the contour 


rotation angle 0. 

In order to solve the integral equations ([2]) for the coupled HcNN — AcNN system we 
consider all conhgurations where the baryon-baryon subsystems are in an S'—wave and the 
third particle is also in an S'—wave with respect to the pair. However, to construct the 
two-body t—matrices that serve as input of the Faddeev equations we considered the full 
interaction including the contribution of the H—waves and of course the coupling between 
the TjcN and AcN subsystems (this is known as the truncated t—matrix approximation 22|). 
This approximation in the case of the NNN system with the NN interaction taken as the 
Reid soft-core potential leads to a triton binding energy which differs less than 1 MeV from 
the exact value 23| . We give in Table [I] the two-body channels that are included in our 
calculation. For a given three-body state (/, J) the number of two-body channels that enter 


/ J 

(*E,,SeJ 

(W,saJ 

(W(Ec),SV(Ec)) (W(A,),'SV(Ac)) 

0 1/2 

(1/2,0),(1/2,1) 

(1/2,0),(1/2,1) 

(1,0) 

(0,1) 

1 1/2 (1/2,0),(3/2,0),(1/2,1),(3/2,1) (1/2,0),(1/2,1) 

(0,1),(1,0) 

(1,0) 

2 1/2 

(3/2,0),(3/2,l) 


(1,0) 


0 3/2 

(1/2,1) 

(1/2,1) 


(0,1) 

1 3/2 

(1/2,1),(3/2,1) 

(1/2,1) 

(0,1) 


2 3/2 

(3/2,1) 





TABLE I: Two-body T,cN channels (^Ecj'SEc); AcN channels (W^SAc)) 3^3^ channels with Sc 
spectator (i7v(Ec))'5 v(Ec))) A'N channels with Ac spectator {iN(Ac)i ^N{Ac)) that contribute to 
a given TicNN — AcNN state with total isospin I and spin J. 
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is determined by the triangle selection rules \ J — \\ < Si < J + \ and \I — ^\ < ii </+|. 
For the parameter b we found that 6 = 3 fm“^ leads to very stable results while for the 
expansion ([T]) we took twelve Legendre polynomials, i.e., 0 < n < 11. 

The two-body interactions are obtained from the chiral constituent quark model of 
Ref. 2^. The NN potentials perfectly describe the S-wave phase shifts 25j] and had also 


been used in the study of the ANN — T^NN coupled channel problem 20, l2l|]. The charmed 


baryon-nucleon potential is derived as explained in Ref. 


26| . For the case of heavy quarks (c 


or b) chiral symmetry is explicitly broken and therefore boson exchanges do not contribute. 
The absence of strange quarks also eliminates the contribution of iF or k exchanges as com¬ 
pared to the strange baryon case. This simplihes a lot the interaction, see Eqs. (8) — (13) 
of Ref. ^ , and gives rise to a potential whose only free parameter would be the harmonic 
oscillator width of the charm quark. We will present our results for different values of be 
to get parameter free predictions. No bound states are found for the charmed two-body 
subsystems. 

There are several facts that should be noted before presenting our results and must be 
considered for our discussion. The kinetic energy associated to the A+ is reduced compared 
with that of the A, what would imply that the A+ would be more strongly bound than 
the A in the case of having identical interactions with the nucleons. However the charmed 


: the strange sector due to the absence of 
8|. Finally, note that he A^ has a positive 


baryon-nucleon interaction is weaker than that o: 
the strange boson exchanges as also noted in Ref. 
charge, whereas the A is neutral. Therefore, Coulomb effects may play a non-negligible role 
in charmed nuclei. In fact, the hypertriton would be unbound if the A were to have a positive 
charge. 

Bearing these considerations in mind, we have hrstly proceeded to solve the Faddeev 
equations for the A^NN and T^^NN systems using the charmed baryon-nucleon and nucleon- 
nucleon interactions derived from the chiral constituent quark model with full inclusion of the 
Ac -f-)- Sc conversion. Let us hrst present the results for the Acd and Sed scattering lengths. 


TABLE II: Acd, ^ 9 ^ 3/2 and Ao^i/ 2 , and Sed, 3^2 1 / 2 ’ scattering lengths, in fm. 


Ao,3/2 

Ao,1/2 

4' 

^ 1 , 3/2 

A' 

1 , 1/2 

-10.27 

9.02 

0.74-Fi0.18 

2.08+iOA7 
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FIG. 2: (a) Fredholm determinant for the J = ll2 and J = 3/2/ = 0 AcNN channels, (b) 
Fredholm determinant for the J = lj2 and J = 3/2 1 = 1 AcNN channels. 


that are shown in Table [TTl Althongh it might be difficult that they can be measured in the 
near future, however they neatly informed us about the possible existence of bound states 
in the different channels. The Scd scattering lengths are complex since the inelastic AcNN 
channels are always open. Both scattering lengths ^i 3 / 2 ’ have a positive real 

part indicating that the interaction is repulsive. The spin 1/2 Acd scattering length, Ao,i/ 2 , 
is also positive. However the spin 1/2 Acd scattering length, ^ 0 , 3 / 2 , is negative, giving rise 
to a bound state with an energy of 271 keV. 

We show in Fig. [2] the Fredholm determinant of the different AcNN (J, J) states ^. The 
1 = 1 channels are repulsive, only the J = 0 channels present attraction. Curiously, as 
already noted in the scattering lengths, the order of / = 0 channels is reversed with respect 
to the strange sector, being the J = 3/2 the most attractive one. This difference can be 

n 

easily understood due to the importance of the A y-)- S conversion in the strange sector [27| . 
The contribution of the A yy S conversion should be even stronger than A yy 77 conversion 
in ordinary nuclei, because it is not suppressed in S'—waves and the A — S mass difference 
is much smaller. When the NA yy AS potential is disconnected, the J = 3/2 channel is 


^ In all / = 0 and 1=1 cases the zero energy corresponds to the binding energy of the deuteron below 
the corresponding threshold, AcNN or HcNN. The deuteron is perfectly reproduced by our model for 
the NN interaction 23. For 1=2 the zero energy corresponds to the HcNN mass because the deuteron 


channel does not contribute (see Table IJ). 
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FIG. 3: (a) Fredholm determinant for the J = 1/2 and J = 3/2 1 = 0 'Ec^N channels, (b) 
Fredholm determinant for the J = lj2 and J = 3/2 1 = 1 EcNN channels, (c) Fredholm 
determinant for the J = 1/2 and J = 3/2 1 = 2 EcNN channels. 

almost not modified while the J = 1/2 losses great part of its attraction. Thus, the ordering 
between the J = 1/2 and J = 3/2 channels is reversed in such a way that the hypertriton 
would not be bound (see Fig. 6(a) of Ref. |20|). The Ac y-)- Sc conversion is less important 
than in the strange sector firstly due to their mass difference, 168 MeV as compared to the 
73 MeV of the strange sector. Besides, it comes reduced with respect to the strange sector 
due to the absence of the strange meson exchanges Sj, giving rise to a smaller NAc ^ VSc 
transition potential. One should also have in mind that in the (/,</) = (0,3/2) channel 
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FIG. 4: Fredholm determinant for the J = \l‘l and J = 3/2 1 = 0 AcNN channels. The solid line 
stands for the full calculation and the dashed line for when the Ac -H- Sc transition is taken to be 
zero. 


the charmed baryon-nucleon interaction with spin-singlet does not contribute, being much 
more repulsive than the spin-triplet one. Note that the ratio for the relative contribution of 
the spin-singlet to the spin-triplet partial waves comes determined from the strange sector 


through the Ap scattering cross section and the hypertriton binding energy 20/ 


We show in Fig. |3] the real part of the Fredholm determinant of the six (J, J) HcNN 
channels that are possible for energies below the Scd threshold. The imaginary part of 
the Fredholm determinant is small. As one can see all channels are repulsive and thus 
uninteresting from the point of view of possible bound states. The larger attraction is found 
in the (/, J) = (1,1/2) channel that in the strange sector presented a quasibound 


state close to the three-body threshold [20| . Such T,NN quasibound state has been recently 
suggested in ^He(K“,7r'^) reactions at 600 MeV/c [28 1. 

To check the relevance of the Ac •H- Sc conversion we have solved the most interesting 
AcNN channels, (J, J) = (0,1/2) and (0, 3/2), switching off the transition between the AcN 
and TjcN subsystems. We plot in Fig. 0] the Fredholm determinant for both cases. The solid 
line indicates the result of the full calculation while the dashed one represents the results 
without Ac O Sc conversion. As can be seen the effect of the T^cNN channel for AcNN is 
not so important as in the strange case. In particular, the order between the two channels 
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is not reversed. As mentioned above, the reduction of the kinetic energy associated to the 
A+ compared with that of the A, implies that the (/, J) = (0,3/2) state is more strongly 
bound than in the strange sector. 

Two hnal remarks are in order. First, there are estimations on the literature about the 
binding energy of heavier charmed hypernuclei {sj and its dependence on the hard core of 
the NAc interaction. Our quark-model approach presents the advantage of having the hard¬ 
core radius hxed by means of quark-antisymmetrization. Second, we have calculated the 
contribution of the Coulomb potential exactly obtaining a contribution of Ec = 131 keV, 
what would give a hnal binding energy of 140 keV for the J = 3/2 charmed hypertriton. 
This result could be easily understood by considering the Anp state as a bound state of a 
deuteron and a A. The Coulomb energy would come given at hrst order by, 

(»l>(r-)|t^(r)|»l>(r-)) 

° {^(f) |S<(r)) ’ ' ’ 

with V{r) = A For small binding energies B, the wave function can be represented by 
^(r) = with k = ■\/2pB, where rj is the reduced mass of the A — d system and B is 
the binding energy. This would give rise to Ec = ak, and using the binding energy B, we 
would obtain Ec = 172 keV, comparable to our exact result and that would not be enough 
as to destroy the bound state. 

To check the dependence of the binding energy of the J = 3/2 charmed hypertriton on the 



FIG. 5: Fredholm determinant for the (7, J) = (0,3/2) AcNN channel for the different values of 
the width parameter for the charmed quark wave function. 
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free parameter of the quark-model charmed baryon-nucleon interaction, we show in Fig. [5] 
the results of the AcNN {I,J) = (0,3/2) state for several models with slightly different 
width parameter for the charmed quark wave function. One should have in mind that in 


Ref. 


291] it was argued that the smaller values of be are preferred to get consistency with 


calculations based on inhnite expansions, as hyperspherical harmonic expansions 


30j, where 


the quark wave function is not postulated. This also agrees with simple harmonic oscillator 
relations he = bn^^- As can be seen the binding energy varies between 1037 and 271 keV, 
and the repulsive Coulomb contribution would vary between 322 and 131 keV, what makes 
the competition between electromagnetic and strong contributions crucial for the existence 
of this state, that would have a binding energy between 715 and 140 keV. 

In summary, we have solved the Faddeev equations for the AeNN and HeNN systems 
using the charmed baryon-nucleon and nucleon-nucleon interactions derived from a chiral 
constituent quark model with full inclusion of the Ac -H- Sc conversion. We present results 
for the binding energy and the Aed and Sed scattering lengths. As compared to the strange 
sector, the kinetic energy is reduced but the interactions are weaker. The smaller contribu¬ 
tion of the Ac O Sc conversion due to the larger mass difference and the smaller transition 
potential reverses the order of the two only attractive channels, (J, J) = (0,1/2) and (0, 3/2), 
the spin-3/2 state becoming the most attractive one. After correcting for Coulomb effects 
the charmed hypertriton would have a binding energy of at least 140 keV. The actual exper¬ 
imental facilities are capable of carrying experiments seeking for theses states what would 
help us in our progress in the knowledge of the baryon-baryon interaction on the heavy-flavor 
sector. 
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